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Abstract: The Blume- Emery-Griffiths model with the dipol-quadrupol 
interaction (£ ) has been simulated using a cellular automaton algorithm 
improved from the Creutz cellular automaton (CCA) on the face centered 
cubic (fee) lattice. The finite-size scaling relations and the power laws of the 
order parameter (M) and the susceptibility (x) are proposed for the dipol- 
quadrupol interaction {£ ). The dipol-quadrupol critical exponent 6i has been 
estimated from the data of the order parameter (M) and the susceptibility 
(x)- The simulations have been done in the interval < i = L/J < 0.01 
for d = D/J = 0, k = K/J = and h = H/J = parameter values on 
a face centered cubic (fee) lattice with periodic boundary conditions. The 
results indicates that the effect of the i parameter is similar to the external 
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magnetic field (h). The critical exponent 5^ are in good agreement with the 
universal value {Sh = 5) of the external magnetic field. 
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The spin-1 Ising model, which is known as the generahzed Blume- Emery- 
Griffiths (BEG) model, can be used to simulate many physical systems. The 
model firstly has been presented for describing phase separation and super- 
fluid ordering in He mixtures [1]. The most general Hamiltonian of the model 
is given by 



Hj = -jJ2SiSj-Kj2SiS] + Lj2iS^S, + S,Sf) + DY,S^ + Hj2Si 

<ij> <ij> <ij> i i 

(1) 

which is equivalent to the lattice gas Hamiltonian under some transforma- 
tions [2] — [4]. (ij) denotes summation over all nearest-neighbor (nn) pairs 
of sites and Si — —1, 0, 1. The parameters J, K, L, D and H are bilinear, 
biquadratic, dipole-quadrupole interaction terms, the single-ion anisotropy 
constant and the external field term. The versions of the model have been 
applied to the physical systems such as the solid- liquid-gas systems [5], the 
multicomponent fluids [6], the microemulsions [7], the semiconductor alloys 
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[8] — [10], He^-He^ mixtures [l,ll],the binary alloys [12] and the magnetic 
spin systems [13] — [16]. 

The spin-1 Ising model with dipole-quadrupole interaction was firstly 
studied using molecular field approximation [6, 17, 18] and the transfer matrix 
method [19] to investigate the tricritical and multicritical points for selected 
values of the I parameter on one and two dimension. The spin-1 Ising model 
has also applied to two dimensional ternary graphite intercalation compounds 
(GIC's) for investigating the influence of the dipole-quadrupole interaction 
on the melting point using Monte Carlo method [20]. The model has been 
used for estimating the critical concentration value of the (GaAs)i_3;Ge2a; 
alloy by Kikuchi approximation [9, 10] and Cluster variation method [21, 22]. 
These studies show that the dipol-quadrupol interaction term is very effective 
on the critical concentration value of the (/// — V)i^xIVx ternary alloys [23] 
and it can be dominant on the phase space [18, 19, 24]. However, the infinite 
lattice critical behaviors of the order parameter (M) and the susceptibility 
(x) with the dipol-quadrupol interaction {£ = L/ J ) has not been studied yet. 
In our previous paper [25] , it was shown that the dipol-quadrupol interaction 
{£) prevents the one of the spin species (5' = -|-1 or —1) on the BEG model 
{J ^ 0, K ^ 0, D ^ 0) and breaks the symmetry in the Hamiltonian 
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similar to the external magnetic field (h) [26]. As a result of this, the dense 
ferromagnetic {df{+), df{—)) and ferromagnetic (F(+), F{—)) phases occur 
on the ground state phase diagram for nonzero £ value. Furthermore, model 
exhibits the reentrant, successive and multi phase transitions in the interval 
—0.15 < £ < 0.15 and it does not show any phase transition at the large £ 
value {£ > 0.15). 

The aim of this study is to investigate the critical behavior of the BEG 
model with the dipol-quadrupol interaction {£ ) and to estimate the power 
laws and the finite size scaling relations for the order parameter (M) and the 
susceptibility (x) in the £ — L/ J > region. In order to expose the effects 
of the dipol-quadrupol interaction {£ — L/J), the biquadratic interaction 
{k — K/J), the single-ion anisotropy term and the external magnetic field 
{h — H/J) have been taken zero { d — 0, k — and h — 0). These terms 
were studied in our previous papers and the effects of them on the phase 
space and the critical behavior of the BEG model were discussed [25] — [29] . 

The temperature dependence of the order parameters (M, Q), the sus- 
ceptibihty (x), the specific heat {C/k) and the Ising energy [Hi) have been 
computed on the fee lattice with linear dimension h— 4, 6, 8, 9 and 12. 
The finite lattice critical temperatures are estimated from the maxima of the 



4 



susceptibility (x) for the foe lattice with periodic boundary conditions. In 
order to expose the influence of the dipole-quadrupole interaction, the ther- 
modynamic quantities have been calculated using CA heating [25], [27] — [29] 
algorithm on the fee lattice for L= 4, 6, 8, 9 and 12 (The total number of sites 
is N — 4L^) with the periodic boundary conditions. The fee lattice can be 
built from four interpenetrating simple cubic lattices. The linear dimension 
L= 12 of face centered cubic lattice corresponds to the L= 20 in simple cubic 
lattice. Hence, the data are analyzed within the framework of the finite-size 
scaling theory for the BEG model with the dipol-quadrupol interaction {£ ). 

The Creutz cellular automaton (CCA) is faster than the conventional 
Monte Carlo method (MC) [31]. The CCA does not need high quality random 
numbers and it is a new and an alternative simulation method for physical 
systems [25] — [29], [30] — [33]. Furthermore, the results obtained using CA 
algorithm and its improved versions are in good agreement with the universal 
critical behavior for the BEG model. Thus the simulations have been carried 
out using a cellular automaton heating algorithm which successfully produces 
the critical behavior of the Ising model [25], [27] — [29]. During the heating 
cycle, energy is added to the spin system through the second variables {Hk) 
after the 2.000.000 cellular automaton steps. The heating rate is equal to 
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O.OSHk per site. The computed values of the thermodynamic quantities 
are averages over the lattice and the number of time steps (2.000.000 ) with 
discard of the first 100.000 time steps during the cellular automaton develops. 
For the finite size lattice, the order parameters M and Q are given by 



M 

L 



^E^^ (2) 



The susceptibility and the specific heat are calculated with 



^^,3(M^ (4) 



C-L (5) 

The expectation values in equation (5) and (6) are averages over the lattice 
and the number of the time steps. 

In order to investigate the critical behavior of the BEG model with the 
dipole-quadrupole interaction, the thermodynamic quantities have been cal- 
culated using CA heating algorithm on the fee lattice for L= 4, 6, 8, 9 and 
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12 (The total number of sites is N — 4L^) with the periodic boundary con- 
ditions. 

For the d — and the k — hne, the model exhibits the second order 
Ferromagnetic- Paramagnetic phase transition [25] for £ = while the sus- 
ceptibility (x) have a characteristic peak. With the increasing £ value in the 
interval < i — L/ J < 0.01, the characteristic peak of susceptibility occurs 
at the higher critical temperature and its peak value decreases. The effect 
of the £ parameter is similar to those coming from the external magnetic 
magnetic field [25, 26, 33]. 

At T — Tc{V), the order parameter values ( M(Tc{V),i) ) have been 
estimated for the increasing i value in the interval 0<£ = L/J<0.01on 
the different lattice sizes L— 4, 6, 8, 9 and 12 and illustrated in figure (la). 
The results show that there is a considerable finite size effect near the i — 0. 
The M{Tc(L),£) tends to zero with the increasing lattice sizes at £ = 0. 
However the M(Tc(L),i) goes to the constant value (~ 0.5) for all lattice 
size with the increasing value of i. The susceptibility data x(Tc(L),£) are 
plotted in figure (lb) for £ values in the interval 0<£<0.01on the different 
lattice sizes L— 4, 6, 8, 9 and 12. Figure (lb) obviously shows that the 
susceptibility data x(Tc(L),£) tend to infinity with increasing lattice size at 
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£ — 0. Around the £ — value, the data show the dependence on the lattice 
size and x(Tc;(L),£) is a decreasing function of £ for the all lattice size. 

The dipol-quadrupol interaction on the spin system shows a similar ef- 
fect with the external magnetic field. As it is seen in figure (1), the order 
parameter (M) and the susceptibihty (x) are the functions of £ . For the 
determination of power-law exponent 5e, the power law relations of the ther- 
modynamic functions ( M and %) for the dipol-quadrupol interaction {£) 
have been considered as similar expressions to the power-laws of external 
magnetic field (h) [26, 34]. Therefore, the power laws of the order parameter 
(M) and the susceptibihty (%) a,t T — Tc have been described by 

M{TcJ) = Be\£\'/'' (6) 

x{Tc,£)^Ce\£f-''^/'' (7) 

where and Ce are the order parameter and the susceptibility amplitudes. 
In the interval < £ < 0.01, the order parameter (M), the critical exponent 
Si{L) and the order parameter amplitude S^(L) on each lattice size (L) are 
obtained from the best fit to straight lines in the Log-Log plot of the data 
(Figure (2a)). The estimated 5^(L) and -B^(L) values are plotted against 
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Lr^/^ in figure (2b) and (c). The extrapolations {L~^l'^ — > 0) of data which 
lie on straight lines give the infinite lattice values as l/5^(oo) = 0.193 ±0.005 
and -B^(oo) = 0.98 ± 0.06. The value 5(, — 5.18 is in good agreement with 
universal value {5h — 5) for external magnetic field [26, 34]. 

In figure 3, the critical exponent 5^(L) and the susceptibility amplitude 
CiiV) are obtained from the susceptibility data (x(7c'(L),£) ). The Log- 
Log plot of x(Tc(L),£) against £ in the interval < £ < 0.001 yields (1 - 
5l(L))/5l(L) and Ci(V) in figure (3a). The extrapolation of these quantities 
gives (1 - 5t{oo))j5t{oo) = -0.80 ± 0.05 and Q(oo) = 0.008 ± 0.001 (Figure 
(3b) and 3(c)). Therefore, the value of ^^(oo) is obtained as 5 which is equal 
to the universal value = 5 for the external magnetic field. 

The finite-size scaling relations [26] of the order parameter (M) and the 
susceptibility (x) related to the dipol-quadrupol interaction I can be defined 
by scaling relations related to the external magnetic field h as 



M = L-^/'^X°(L'5^^/^ |£| , L^/'^e) 



(8) 



(9) 



with £ = (T - Tc{oo))/Tc{oo) [26, 34]. 
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At T — Tc{oo), the scaling functions X° and Y° are asymptotically re- 
produced as 



where x =]J^PI'^ and are the order parameter and the suscepti- 

bility amplitudes. The finite-size scaling plots of the order parameter data 
M{Tc{oo),€) at the infinite lattice critical temperature {e — Q ) are illus- 
trated in figure 4. For (5^ = 5, /3 = 0.31 and v — 0.64 universal values, the 
order parameter data lie on a single curve with the slope equal to 1/5^ = 0.2 
in the interval 0.229 < a; < 2.875. The estimated value is in good agree- 
ment with the value of the field critical exponent ( 5h—b). Furthermore, the 
order parameter amplitude in equation (9) is estimated as Bg^ = 1.096. This 
value is in agreement with the hnear extrapolation result (Figure. 2(c)). The 
finite-size scaling plots of the susceptibility data x{Tc{oo),€) are shown using 
(5^ = 5, 7 = 1.25 and v = 0.64 universal values in figure 5. As it is seen in 
figure 5, the scaling susceptibility data in the interval 0.236 < x < 3.018 lie 
on a straight line with the slope (1 — Si)/Si = —0.8. It gives the value of the 
dipol-quadrupol critical exponent 5^ as 5. The value of the straight line at 



X°(x) = Bex^/^^ 



(10) 



Y°{x) = Cix^^-^'^l^' 



(11) 
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x — gives the susceptibility amplitude as Q = 0.063. 

The previous studies implied that the dipol-quadrupol interaction i on the 
spin-1 Ising model Hamiltonian is considered as a magnetic field-like pertur- 
bation [8, 35]. Our simulations exposed that the dipol-quadrupol interaction 
parameter i has prevented the one of the spin species {S = +1 or —1) and it 
has broken the symmetry in the Hamiltonian. Really, the £ acts similar to the 
external magnetic field h over the spin system. However, the order parameter 
M{Tc, i) and the susceptibility x(Tc, i) behave according to the equations of 
state ( M ~ i^/^ and x ~ e^'^'^'^/^^) in the interval < £ = L/J < 0.01 at the 
critical temperature. Moreover, the order parameter (M) and the suscepti- 
bility (x) data are scaled well within the framework of the finite size scaling 
hypothesis (at Tc{oo) obtained for i — 0). The obtained value of the critical 
exponent 5^ are in good agreement with the universal value {5h — 5) of the 
external magnetic field. These results indicate that the dipol-quadrupol in- 
teraction parameter i is similar to the external magnetic field h over the spin 
system. Therefore, i can be considered as a magnetic field-like perturbation. 
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Figure Captions 

Figure. 1. At T = Tc(L) on the different lattice sizes L= 4, 6, 8, 9 and 
12 in the interval < £ < 0.001 for d = and k — (a) the order parameter 
values M against i, the susceptibility values x against £. 

Figure. 2. On the difTerent lattice sizes L— 4, 6, 8, 9 and 12 (a) Log-log 
plots of M versus £, (b) The plot of 5e{L) versus L'^^", (c) The plot of Be{L) 
versus L~^/'^ . 

Figure. 3. On the different lattice sizes L— 4, 6, 8, 9 and 12 (a) Log-log 
plots of X versus £, (b) The plot of ^^(L) versus L~^^'^, (c) The plot of C^(L) 
versus L"^/*^ . 

Figure. 4. The finite-size scaling plot of the order parameter (M) for 
T < Tc{oo) on the different lattice sizes L= 4, 6, 8, 9 and 12 
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Figure 5. The finite-size scaling plot of the susceptibihty (x) for T < 
Tc{oo) on the different lattice sizes L= 4, 6, 8, 9 and 12 
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